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For a class of generalized Hubbard models, we determine the maximal stability region for the 
superconducting ry-pairing ground state. We exploit the Optimized Ground State (OGS) approach 
and the Lanczos diagonalization procedure to derive a sequence of improved bounds. We show that 
some pieces of the stability boundary are asymptotic, namely independent on the OGS cluster size. 
In this way, necessary and sufficient conditions are obtained to realize superconductivity in terms 
of an 77-pairing ground state. The phenomenon is explained by studying the properties of certain 
exact eigenstates of the OGS hamiltonians. 



Generalized Hubbard models are important theoreti- 
cal frameworks for the study of superconductivity. Apart 
from special cases, they are not solvable and rigorous re- 
sults on their physical properties are quite valuable. 

The inclusion of nearest neighbour extra couplings be- 
yond the on site Coulomb interaction has a long history 
and is still an interesting issue. Indeed, the qualitative 
effects of these interactions is not definitely understood 
and examples can be provided where they are not negli- 
gible at all. 

As is well known, a good marker for superconductivity 
is off-diagonal long-range order (ODLRO) jl|, a property 
that makes sense in any number of dimensions and im- 
plies both Meissner effect and fiux quantization. Ground 
state ODLRO can be detected by studying the asymp- 
totic behaviour of fermion correlation functions ||^. Of 
course, if the ground state is analitically known, it can 
be checked explicitely. This is the case of the so-called 
77-pairing ||] states that exhibit ODLRO and, under some 
constraints, can be the ground states of certain general- 
ized Hubbard models. 

When an 77-pairing state is discovered to be an ex- 
act eigenstate, the next problem is to determine the re- 
gion in the coupling space where it is also the ground 
state. To answer this question many analytical meth- 
ods have been developed to establish rigorous bounds 
for the superconducting region. Among them, we re- 
call the positive semidefinite operator approach and 
the bounds derived by application of Gerschgorin's theo- 
rem [0 . The algorithm which however appears to be the 
simplest and most powerful is the Optimal Ground State 
(OGS) scheme proposed for generalized Hubbard mod- 
els and recently applied to the case of next to nearest 
neighbour couplings |^ . The method is based on the ex- 
act diagonalization of a certain local hamiltonian defined 
over a cluster of sites. If the cluster is made larger, the 



superconducting region is generally expected to expand. 
In the limit of an infinite cluster we obtain exact bounds. 

For simplicity, in the following we shall call supercon- 
ducting (SC) region, the subset of coupling space where 
the r^-pairing state is the ground state. In this Letter, 
we apply the OGS algorithm to study the stability of the 
superconducting ry-pairing state with momentum tt. We 
discuss the behaviour of the OGS bounds as a function 
of the cluster size using the Lanczos algorithm to diag- 
onalize the cluster hamiltonian. We obtain an improved 
SC region that can be considered numerically asymptotic 
and discuss in details the inclusion problem by stating 
the conditions under which larger clusters are expected 
to give better bounds. Another interesting result is that 
some pieces of the boundary between the SC and non 
SC regions are independent on the cluster size and de- 
termine necessary and sufficient conditions for 77-pairing 
superconductivity. We explain these stable boundaries 
by means of certain exact eigenstates of the OGS hamil- 
tonians whose properties are crucial in this respect. 

Let us consider the Hamiltonian of a one dimensional 
generalized Hubbard model (we denote by the sum 
over nearest neighbour sites) 

H = -t {c\aC3<y + c]„C^a) + 

{^,3),'T=-\,i 

+ X ^ (7^j,_^ + 77j,_^)(ct^Cj„ -I- c]^C,^) -f 

+ y ^ (77, - 1)(7., -1)+YY, (pIp, + p]p,) + 
{i.j) (ij.) 

~^2'^xy{SfSj + Sj' ) + JzSf Sf , (1) 
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where Cia and cj^. are canonical Fermi operators obey- 
ing {cL>Cja'} = SijSaa- and {cia,Cj„'} = {cl^,c]^,} = 0. 
The number operators are defined as usual: riio- = Ci^- 



and Tii 



nil- The operator creates pairs 



= c]^|cj^. The Heisenberg exchange interaction is writ- 
ten as usual in terms of su{2) operators — cl^Cii, 
S- = (5'+)t and = i(n| - n^). 

The Hamiltonian in (m contains many couplings: X 
parametrizes the bond-charge repulsion interaction which 
has been related to high-Tc materials ; U is the usual 
on site Coulomb interaction; V is the nearest neighbour 
charge-charge coupling; Y controls the pair hopping term 
as in the Penson-Kolb- Hubbard models . Finally, Jxy 
and Jz are the t— J like Heisenberg exchange couplings. 
We introduce the 77-pairing operator 



from which we build the state 



|,y) = (^t)A'/2|o), 



(2) 



(3) 



where |0) is the empty state and N is the number of lat- 
tice sites. The state \ri) is an eigenstate of H provided 
V + 2Y = 0. In this case it describes a half-filled state 
with energy 



+ ^-iU + AV), 



(4) 



and can be shown to possess ODLRO. Since £'+ is an up- 
per bound for the ground state energy, a strategy to proof 
that 1 77) is the ground state is to find a lower bound 
and a region in the coupling space where E- The 
same procedure applies also to other exact eigenstates 
like, for instance, the 77-pairing state with zero momen- 
tum. Also, ii V ~ Y = Jxy ~ Jz ^ and t ~ X, exact 
eigenstates can be built by applying 77^ to eigenstates of 
the U = +00 standard Hubbard model due to the fact 
that t ^ X implies a conserved number of doubly oc- 
cupied sites. However, we remark that the state [77) is 
particularly interesting from a methodological point of 
view as can be seen by the complete failure of the naive 
Gerschgorin approach jl^ . 

Lower bounds for the ground state of H may be ob- 
tained following the OGS approach The Hamilto- 
nian (|l|) is written as 



H 



00 

E 



i (2) 
''n,n+ 



1), 



(5) 



where h 

.(2) 



(1) 



contains operators acting only on site n and 
^rTn+i links site n to site 71, -I- 1 and depends on operators 
acting on both. To recast (H), we introduce extended 
operators 



n+k-2 



m—n 

n+k-2 



(6) 



m—n+l 



for any integer k > 2. The local hamiltonian /i'^'^^ de- 
scribes a cluster of k sites. Like H, also /i^*^^ is symmetric 
under 77'''. All the states 



k sites 



p integer 



(7) 



are degenerate with energy E^ and are precisely those 
needed to build the \ri) state on the infinite lattice (see ||] 
for a complete discussion of the k = 2 case). The Hamil- 
tonian can be written in terms of ft,*^*^) as 



H= y: hi''- 



(8) 



If we denote by Eq{N) the ground state energy for a sys- 
tem of N sites, the asymptotic ground energy per site is 
by definition 



£q = lim 



EojN) 
N '■ 



and, for each fc, satisfies the rigorous bound 



1 



k-l 



(9) 



(10) 



where (j{A) denotes the spectrum of the operator A. The 
normalization factor l/(fc — 1) takes into account the 
number of terms in Eq. (^) which contain a given site. 
The right hand side of Eq. ( |l0|) depends on the cluster 
size k and a better bound is expected as k increases. 
However, strictly speaking, this is false. Let us write a 
cluster of k sites in terms of smaller clusters 



From (O) we obtain the exact inequalities 



> 



k-l-l 



1 



[k-i) 



I 



1 



and in particular 



(2fe-l) 



(k) 



(11) 



(12) 



(13) 



This allows to build sequences of bounds converging to 
the exact bound in the infinite cluster size limit. In more 
details, Eq. ( p^ splits into disjoint sequences of cluster 
sizes as follows 



2c3c5c 
4c 7c 13c 



(14) 
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where the notation means that each sequence gives bet- 
ter and better bounds. We notice that in general two 
sequences are not related at all and, in particular, 3 C 4 
may be false as we shall see in explicit examples. What 
can be stated in full generality is that the minimal choice 
fc = 2 is always the worst bound since from 



(JV+l) 



> 



A^- 1 



N 



'0 



1 



'(2) 



we proof inductively that VA'^ > 2 we have 



(2) 



(15) 



(16) 



Keeping these remarks in mind, we study the size depen- 
dence of the conditions under which (|^) with Y = —2V 
admits {r]) as its ground state by explicit diagonalization 
of h^^^ on clusters of increasing sizes. The OGS method 
requires diagonalization of the local hamiltonian in all 
sectors of definite up and down electron numbers; for the 
numerical diagonalization we use the Lanczos algorithm. 
In the following we shall always assume t = I and denote 
by L the cluster size. 

Since the couplings constants space is large, we decide 
to discuss separately what happens with the Heisenberg 
exchange interaction switched on or off. Let us begin 

with Jxy = Jz =0. 

In Fig. 1 we show dX X — the size dependence of 
the bounds when V > —1. As can be seen, there are 
regions where the corrections are definitely negligible be- 
yond L = 3, i.e. V > —0.4. On the other hand, around 
F = — 0.5, size effects can be important up to large clus- 
ter sizes. We remark that this figure does not show any 
non trivial relationship among the bounds obtained at 
different L: they just improve monotonically. 

In Fig. 2 we show the best bounds obtained with L = 6 
at several values of X. In the inset we expand the region 
around V = —2. A remarkable feature of the plot is that 
an enveloping straight line appears around V — ~1. Suc- 
cessive corrections are quite small and the region shown 
can be considered maximal from any practical point of 
view []l3| . 

In Fig. 3 we plot at four different X the difference 
AU{L) = U{L) — U{2) between the boundary curves at 
L > 2 and the minimal one at L — 2. The inclusion 
tree (|lj) is non trivially satisfied and indeed the L — A 
bound is not always better than the L — 3 one. As a sec- 
ond remark, we observe that at each X there is a piece of 
the boundary where finite size corrections vanish. This 
turns out to happen between two of the L — 2 boundary 
points. As shown in the result at L = 2 is that a 
sufficient condition for to be the the ground state is 



y < 0, 



U <-2 max { 2 + 2V, 2|1 - 2X\ + 2V,V - 



(17) 



V 



For < X < 1 (we study this case only), the differ- 
ence AU vanishes between the intersections of the curves 
U = -4(1 + V) and U ^ ~2{V ~ [t - X f/V), namely 
for 1| < ^X{2-X). 

Let us now discuss why this stable boundary subset ap- 
pears. For each value of L, the normalized cluster Hamil- 
tonian ■j^zih'^^^ has many eigenstates \e['^\u,V, X)) 
(i = 1, . . . , dim(/i)) which we label by their eigenvalue. 

Let X play the role of a parameter; following the OGS 
approach, the inequalities -Ej-^'' > i?+ determine the su- 
perconducting region in the (L/, V) plane. Each point 
of its boundary satisfies E^'"^ = E+ for some index i. 
Hence, if a subset of the boundary turns out to be L 
independent, a possible reason can be the existence of 
an eigenvalue independent on L. A trivial case is pro- 
vided by the states (?7^)p|0) {p integer) where |0) is the 
empty state for h^'"^ However, in this case, the condi- 
tion = E+ is identically satisfied for all J7, V and 
X and does not determine any boundary. To find a non 
trivial eigenstate with eigenvalue independent on L we 
can consider the one particle sector (i.e. n| = 1, = 
or viceversa). The two states 



|0), 



(18) 



are indeed exact eigenstates of -j^zxh!^^^ provided U — 
—4(1 + V) and in this case their eigenvalue is precisely 
E+ = — 1 (i = 1). The states l^o-) are thus responsible 
for the stable boundary. To understand why it is con- 
fined to |y -I- 1| < \/2X — we introduce additional 
eigenstates of Indeed, on the line U = —A{V + 1), 

the su{2) singlet state {X ^ 1) 



l7) 



Si 



l-(-l) 



i+L 



-A ) |0), (19) 



can be shown to be an exact eigenstates of with 
eigenvalue E+ if and only if p = (2 -f y)/(l — X) and 
V ^ -l± V2X - X2. This is the L > 2 generalization 
of the state \'4'±) discussed in [|j in the L = 2 case. It 
forbids to extend the bounds associated to \Sa) beyond 
the points \V + l\ = V2X - X^. The case X ^ 1 is 
singular and must be treated separately; the number of 
doubly occupied sites is conserved and splitting may oc- 
cur. For instance, the state I7+) with V = splits into 
the independent eigenstates |7i) =pl\0). 

The above analytical and numerical results lead us to 
the conclusion that the inequality U < —4(1-1-1^) is a nec- 
essary and sufficient condition for \r]) being the ground 
state in the subset of the coupling space constrained by 
the conditi ons t = 1, Y + 2V = 0, 0<X<1 and 
\V + 1\ < ^X{2-X). 
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The above facts do not change quahtatively when the 
Heisenberg exchange interaction is switched on. To sim- 
pHfy the analysis, we consider the special point Jxy = 
Jz = —2Y which allows for a comparison with the re- 
sults illustrated in . 

The plot of Fig. 4 is analogous to that of Fig. 1; we 
show the OGS bounds in the {U, V) plane a,t X = 0. A 
part of the boundary is clearly independent on L. This 
fact is clearly visible in the left inset and mostly in Fig. 5 
where (as in Fig. 3) the asymptotic part of the bound- 
ary can be seen as a function of X. In particular, the 
left edge is at y = — 1 and is independent on X. These 
features can be analyzed as in the previous case. The 
bounds valid at L = 2 are 



V<0, 

U<-2 max! 0, 2 + 2F, 2|1 - 2X| + 2V, 4V 



(20) 

V 



and the stable part (with < X < 1) is that between the 
intersections of the lines U = 0, U = —4(1 + V) and the 
cuTveU = -2{4V-{1-Xf/V). Thus wc obtain the fo l- 
lowing interval for V: -1 <V < ^{1- V3-AX + 2X^). 
As before, the Hne U = —4(1 + V) appears at all L since 
on it the states \Sa) are eigenstates of j^h^^^ ■ More- 
over, as before, there are two states which forbid to cross 
the above interval for V. At the right edge, such a state 
has the form of Eg. (pj|) wit h p = 2(1 - V)/{1 - X) and 
y = i(l - V3-4V+ 2A:2). At the left edge, V = -1, 
the exact {X independent) eigenstate is instead 



IV) 



i>3 



|0). 



(21) 



To summarize, in this case, our conclusion is that un- 
der the conditions t=l, Y + 2V = 0, 0<X<1 

^(1 



and -1 < V < 



2^-1- -\/3 — 4:X + 2X^), a necessary 
and sufficient condition for \ri) being the ground state 
is f7 < -4(1 + y). 

To conclude, in this Letter we have considered gener- 
alized Hubbard models with nearest neighbour couplings 
and the problem of determining when the ground state is 
a superconducting 77-pairing state. By diagonalizing local 
hamiltonians associated to clusters of sites with different 
sizes we studied the convergence of the OGS bounds. As 
predicted in 1^ , it may happen that the bounds obtained 
with the smallest clusters are actually exact. This pe- 
culiar situation seems rather typical and indeed we have 
shown that there exist subsets of the bounding region 
which are asymptotic and remain unchanged as the clus- 
ter size is varied. We clarified the origin of the phe- 



nomenon by providing several exact eigenstates which 
play a crucial role in its derivation. 

We thank A. Montorsi for very useful discussions on 
generalized Hubbard models and 77-pairing superconduc- 
tivity. 
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FIG. 2. Best OGS bounds obtained with L = 6 clusters. The Heisenberg interactions are switched off Jxy = Jz = 0- The 
different curves correspond (from bottom to top) to X = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 1.0. 
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FIG. 3. Existence of an asymptotic boundary. The plots show U{V; L) — U{V; 2) as a function of V (always in units of t) 
for four values of X. The function U{V; L) is the curve obtained from the OGS bounds using clusters of L sites. The inset at 
X = shows a non trivial inclusion tree as the cluster size is increased. See also Fig. 4. 
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FIG. 5. As Fig. 4, but with J^y = = —2Y. The inset magnifies a portion of the curves to confirm that the OGS bound 
not asymptotic for F < — 1. 



